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Abstract. We use improved Monte-Carlo algorithms to study the antiferromagnetic 2D-Ising 
model with competing interactions Ji on nearest neighbour and J2 on next-nearest neighbour 
bonds. The finite-temperature phase diagram is divided by a critical point at J2 — Ji/2 where 
the groundstate is highly degenerate. To analyse the phase boundaries we look at the specific 
heat and the energy distribution for various ratios of J2 / Ji . We find a first order transition for 
small J2 > Ji/2 and the transition temperature suppressed to Tc = at the critical point. 



1. Introduction 

The antiferromagnetic J1-J2 spin-1/2 Heisenberg model [1-4] and related bosonic models on the 
square lattice [5-8] with hopping ti and repulsion terms Ji on nearest neighbour (NN) and next- 
nearest neighbour (NNN) bonds have been intensively studied over recent years via different 
approaches. The treatment of the bosonic model via Quantum Monte-Carlo (QMC) algorithms 
suffers from freezing problems in the intermediate regime at J2 = Ji/2 which can be traced to 
a groundstate degeneracy of the static limit with hopping terms = 0. This motivated us to 
restudy the 2D-Ising model with competing interactions and solve these freezing problems via 
improved classical Monte-Carlo (MC) algorithms. The 2D-Ising model has a long history and is 
still a topic of discussion [9-14]; especially the character of the phase boundary for J2 > Ji/2 
is an open question [15-18]. There are two scenarios discussed: Landau and Binder as well as 
Malakis et al [14,16] used MC methods to infer a second order transition for all J2 > Ji/2 with 
non universal exponents. In contrast Lopez et al and Anjos et al use variational methods [15, 18] 
and find a first order transition for small Ji ^ J2 > >^i/2 and a continuous transition only for 
larger J2. Our simulation strengthens the latter scenario [19]. 

We outline the model and the methods in section [5] and O and present our results in section 
m We end with a short discussion and outlook. 

2. Model 

The model of hardcore bosons on a square lattice with competing interactions on NN bonds and 
NNN bonds is described by the hamiltonian: 

-f^boson = -h ^{b\bj + b]hi) - t2 {b\bj + b]bi) + Vi ^ n^n^- + V2 ^ n^n^- (1) 
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Figure 1. Sketched are both ordered phases 
and a third groundstate configuration at 
J2 = «/i/2 (middle). Shaded areas mark 
flipped hnes. 
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Figure 2. Binder cumulants depending on 
temperature and L for J2 = Ji- In the inset 
the intersection area of the cumulants with 
errorbars gives Tq/Ji = 2.082 ± 0.005. 



This hamiltonian can be mapped onto an eqivalent spin model which conincides with the J1-J2 
spin-1/2 anisotropic Heisenberg model apart from the sign in front of the hopping terms tj. 
Taking the static limit = the hamiltonian reduces to that of the Ising model with competing 
antiferromagnetic interactions Jj > 0: 

-falsing = Jl ^ SiSj + J2 ^ SiSj , Si = ±1. (2) 
NN NNN 

We will study square lattices of linear extent L with periodic boundary conditions. 
For J2 < Ji/2 the groundstate of ^ is Neel-ordered (see Fig. [U left) with energy: 

^Neel = -2iV( Jl - J2) , N = L\ (3) 

For J2 > Ji/2 the system orders in the collinear (or superantiferromagnetic) phase (see Fig. [H 
right). In this state every spin has two parallel and two antiparallel aligned nearest neighbours. 
Thus, the energy depends only on the antiparallel diagonal bonds: 

Ecou = -2NJ2 . (4) 

At the critical point J2 = Ji/2 separating these two phases the transition temperature is 
suppressed and the groundstate is highly degenerate of order 2^+^ — 2 (one of these states 
is shown in the middle panel of Fig. [T]). Due to this groundstate degeneracy the free energy 
is characterised by many local minima. This leads to freezing problems for simple importance 
sampling MC simulations. 



3. Methods 

In order to compute the transition temperature Tc for ratios of J2/J1 very close to the critical 
point J2 = Ji/2 we used parallel tempering [20-23] MC methods with additional line flip 
updates [19]. 

To detect the flnite temperature phase transitions we use the Binder cumulant (reduced 
fourth order cumulant) U4 = 1 — ^^^i^i [14, 24, 25] of the order parameter M{q) = S{q)/N , 
where 

S{q) = l^Y.^''-^''~''^{S,S,) (5) 
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Figure 3. Plotted is Tq over the ratio 
J2/J1. Data were produced by parallel 
tempering MC simulations. 
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Figure 4. Energy histogram for several 
ratios of J2I J\ at Tq depending on given L. 



is the respective structure factor (for an example see Fig. [2]). 

To determine the critical exponent v we calculate the derivative of C/4 with respect to T which 
is proportional to L^l^ [14]. 

The behaviour of the specific heat and energy histograms [26-29] for various ratios of J2I J\ 
provides further details to analyse the character of the given finite temperature phase boundaries 
and the transition at the critical point J2 = Ji/2. 

4. Results 

In Fig. [3] we present the finite temperature phase diagram for the classical Ising model (solid 
line) for various ratios \ J2l J\ — 0.5| > 0.005 very close to the critical point. The data is in good 
agreement with MC results from Landau and Binder for | J2/J1— 0.5| > 0.1 [10,12] and for J2 = J\ 
{Tc/Ji = 2.082 ± 0.005 as shown in Fig. ED with the recent result Tc / Ji = 2.0823 ± 0.0017 [16]. 

Calculating lattice sizes up to L = 300 we find the critical exponent v = 1.002 it 0.010 for 
J2 = O.3J1. This and a slowly emerging peak in the specific heat confirms the assumption that 
for J2 < Ji/2 the phase transition lies in the Ising universality class [30,31]. For J2 > Ji/2 
the character of the phase transition is an open subject. We used a histogram method [26-28] 
to verify a first order transition [29] for small J2 > Ji/2. The slowly emerging double peak 
in the energy distribution for J2 = 0.6 Ji and J2 = 0.65 Ji in Fig. H] compared to the single 
gaussian peak for J2 < Ji/2 indicates a first order transition. For larger J2 the shape of the 
energy distribution only shows a broadened peak for the simulated lattice sizes. The calculation 
of V for J2 > Ji/2 strongly depends on the lattice size because of large crossover scales near 
the critical point and above. For J2 > Ji we find the critical exponent u k, 1 {v = 1.02 it 0.02 
for J2 = 1.2 Ji, V = 0.96 it 0.02 for J2 = 1.5 Ji) for large lattice sizes which indicates that 
an Ising-like second order phase transition is recovered for J2 > Ji. So we believe that the 
character of the phase transition varies with J2 and is transfered from a first order transition 
for Ji/2 < J2 < Ji to a continous transition for larger J2 with universal exponents. 

At the critical point J2 = Ji/2 the finite size scaling analysis of the specific heat [19, 27] 
suggests a suppression of the transition temperature to Tq = 0. 

5. Discussion 

We used substantially enhanced MC methods to determine transition temperatures and 
characterise the phase boundaries of the frustrated 2D-Ising model. We solved the freezing 



problems which arise near the critical point at J2 = Ji/2 due to a groundstate degeneracy 
and calculated critical temperatures in the direct vicinity. We found the phase transition for 
<^2 < to be continous and Ising-like. For J2 = Ji/2 the phase transition is suppressed to 
zero temperature. The character of the phase boundary on the right hand side of the phase 
diagram was discussed via histograms. We decide the discussion about the character in favour 
of a first order transition for small J2 > Ji/2 and find strong indications that the transition is 
transfered into a continous one with universal Ising exponents for large J2 > Ji- However, this 
is complicated to show via MC simulations because of large crossover scales and the need for 
very large lattices. 

Next we want to introduce finite hopping terms ti and study the bosonic model given in ([T]) 
using parallel tempering QMC simulations [32]. 
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